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COMPUTATIONAL TRANSITION AT THE UNIQUENESS THRESHOLD 

ALLAN SLY 

Abstract. The hardcore model is a model of lattice gas systems which has received much atten- 
tion in statistical physics, probability theory and theoretical computer science. It is the probability 
distribution over independent sets I of a graph weighted proportionally to A'^' with fugacity pa- 
rameter A. We prove that at the uniqueness threshold of the hardcore model on the d-regular 
tree, approximating the partition function becomes computationally hard on graphs of maximum 
degree d. 

Specifically, we show that unless NP=RP there is no polynomial time approximation scheme for 
the partition function (the sum of such weighted independent sets) on graphs of maximum degree 
d for fugacity Ac(d) < A < Ac(d) + e{d) where 

^ ^{d-iy-' 



(d-2)'* 

\^ ' is the uniqueness threshold on the d-regular tree and e(d) > is a positive constant. Weitz [34) 

\^ ■ produced an FPTAS for approximating the partition function when < A < Ac(d) so this result 

jy^ ' demonstrates that the computational threshold exactly coincides with the statistical physics phase 

O ■ transition thus confirming the main conjecture of [28]. We further analyze the special case of 

A = l,d = 6 and show there is no polynomial time approximation scheme for approximately 

counting independent sets on graphs of maximum degree d — 6, which is optimal, improving the 

^ I previous bound of d = 24. 

Our proof is based on specially constructed random bi-partite graphs which act as gadgets in 

a reduction to MAX-CUT. Building on the involved second moment method analysis of [28] and 

|/-\ ■ combined with an analysis of the reconstruction problem on the tree our proof establishes a strong 

■ I version of "replica" method heuristics developed by theoretical physicists. The result establishes 

^— s ■ the first rigorous correspondence between the hardness of approximate counting and sampling with 

(^ I statistical physics phase transitions. 
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1. Introduction 

The hardcore model is a model from statistical physics representing hardcore interaction of 
gas particles. It is a probability distribution on independent sets I of a graph weighted as ^X'^' 
where A is a positive parameter called the fugacity and Z is a normalizing constant called the 
partition function. Physicists and probabilists have done extensive work towards identifying the 
phase transitions and other properties of the model. 

In computational complexity approximately counting (weighted) independent sets is a central 
problem. The hardcore model is of key importance as this is exactly the problem of producing 
an FPRAS (fully polynomial randomized approximation scheme) for Z, the partition function. 
When A is small the hardcore model has rapid decay of correlations and the partition function 
can be approximated either using MCMC or through computational tree methods [34]. For larger 
fugacities long range dependencies may appear and the problem is known to be hard when A is 
sufficiently large. 

In this paper we determine a computational threshold where approximating Z becomes hard. 
Using an ingenious computational tree approach Weitz [34] produced a PTAS for approximating Z 
when A < Xc{d) where 

A (d) - (^-^)'"' 

^'=^^^- {d-2r 

is the uniqueness threshold for the hardcore model on the infinite d-regular tree [TTj, the point 
at which long range dependencies become possible. Mossel, Weitz and Wormald [28] showed that 
beyond this phase transition local MCMC algorithms fail and conjectured that it gives the threshold 
for computations hardness. While such statistical physics phase transitions are believed to coincide 
with the transition in computational hardness of approximating the partition function for a number 
of important models no such examples had been proven. Our main result essentially confirms the 
conjecture of |28j giving the first such rigorous example. 

Theorem 1. For every d > 3 there exists e{d) > such that when Xc{d) < A < Xdd) + e(d), 
unless NP=RP, there does not exist an FPRAS for the partition function of the hardcore model 
with fugacity X for graphs of maximum degree at most d. 

While we believe the result holds for all A > Ac, for technical reasons (specifically showing that an 
explicit function of three variables attains its maximum at a prescribed location, see Section 11.3.11 
for details) the result is limited to A close to criticality. This limitation notwithstanding, it clearly 
demonstrates the central role played by the uniqueness threshold. 

When A = 1 the hardcore model is simply the uniform distribution over independent sets and the 
partition function is simply the number of independent sets and as such this case is of particular 
interest. When d < 5 Weitz's result provides a FPRAS as Xc{d) > 1. Conversely it is known 
that with d > 25 the problem is computationally hard [9]. While the case d = 6,A = 1 does not 
fall within the scope of Theorem [U using a computer assisted proof, we establish the necessary 
technical condition and prove the following result. 

Theorem 2. Unless NP=RP for every d > 6 there does not exist a fully polynomial approximation 
scheme for counting independent sets on graphs of maximum degree at most d. 

1.1. Background and Previous Results. Even on graphs of maximum degree 3 the problem 
of exactly counting independent sets is #P hard [13] and as such one can at most ask when it is 
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possible to approximately count independent sets, that is when an FPRAS exists. As the model 
is self-reducible, approximate counting is equivalent to approximately sampling from the partition 
function [29]. This has led to a major line of research in analyzing the performance of MCMC 
techniques, particularly the Glauber dynamics. 

When A < ^^^ the Glauber dynamics mixes rapidly [20j which in particular gives an FPRAS for 
counting independent sets on graphs of maximum degree at most 4 (see |10j for similar bounds). 
Weitz [34j showed that the hardcore model has a decay of correlation property called strong spatial 
mixing whenever A < Ac which implies rapid mixing on graphs of sub-exponential growth. Moreover, 
his paper gives a deterministic polynomial time approximation scheme on all graphs when A < Ac 
through a computational tree approximation. 

Finding the ground state of the hardcore model, the largest independent set, is of course a 
canonical NP-hard problem and is hard to approximate even on regular graphs of degree 3 [3]. 
Intuitively the problem of counting becomes harder as A grows as this places more mass on the 
larger, harder to find, independent sets and indeed such hardness results have been established. 
In [20] it was shown that there is no FPRAS (assuming NP=RP) when A < c/d for c ~ 10000. 
In the case of A = 1 this was improved to d > 25 in [9] using random regular bi-partite graphs as 
basic gadgets in a hardness reduction. They further showed that with high probability the mixing 
time of the Glauber dynamics on a random bipartite d-regular graph is exponential in the size of 
the graph. Calculations of [9] led the authors there to speculate that Ac may be the threshold for 
hardness but the evidence was not conclusive enough to make such a conjecture. 

1.1.1. Replica Heuristics. The replica and cavity methods and heuristics have provided powerful 
tools (often non-rigorous) in the study of a wide range of random optimization problems and 
predictions for the behavior of spin glasses and dilute mean fields spin systems [22l[23]. Developed 
by theoretical physcicits, in in some cases these heuristics have been made rigorous, notably the 
SK model [33], solution space of solutions to random constraint satisfaction problems fl] and the 
assignment problem [2]. In dilute spin glass models such methods have given rise to powerful new 
algorithms such as survey propagation (see e.g. [T9]). 

Random regular bi-partite graphs are widely known to be locally tree-like with only a small 
number of short cycles. The statistical physics theory makes the following predictions for the 
hardcore model on typical random bi-partite d-regulars. The first is that the model is expected to 
exhibit spontaneous symmetry breaking for A > Ac. When A < Ac correlations decay exponentially 
and the configuration (independent set) is essentially balanced between the two halves of the bi- 
partite graph. By contrast when A > Ac the configuration separates its mass unevenly placing VL{n) 
more mass on one side or the other. Configurations with a roughly equal proportion of sites on each 
side make up only an exponentially small fraction of the distribution. This is intuitively plausible 
as the largest bi-partite sets will be those containing most of one side of the graph or the other. 

The second is that this symmetry breaking splits the configuration space into two "pure states" of 
roughly equal probability. We will denote the "phase" of the configuration as the side of the graph 
with more sites. Conditional on the phase the spins of randomly chosen vertices are assumed to be 
asymptotically independent and the local neighbourhood of the configuration are given by extremal 
measures. This conditional independence is a crucial element of cavity-method type arguments. 

A first moment analysis of [9] suggested that configurations obey the first prediction but their 
proof proceeded without specifically proving it. In a technical tour de force the prediction was 
rigorously established for Xc{d) < A < Ac(d) + £{d) in [28j using an involved second moment 
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method analysis together with the small graph conditioning method. The restriction to the region 
A < Xc{d) + e{d) is somewhat surprising at first as the problem ought to become easier as A grows. 
It is the result of a technical difficulty in estimating the second moment bound. Even establishing 
this for A close to the critical value took up fully a third of the proof. As a central part of our proof 
is a modification of this method the same restriction applies. 

Based on establishing the symmetry breaking [28] showed that any local reversible Markov Chain 
has mixing time exponential in the number of vertices by establishing a bottleneck in the mixing 
on asymptotically almost all random d-regular bi-partite graphs. This bound is tight as subsequent 
results [27^ Theorem 4] imply rapid mixing on almost all random bi-partite graphs when A < Xc{d). 
Based on these finding they made the following conjecture. 

Conjecture 1.1. f j28j j Unless NP=RP for every d > 4 and Xc{d) < A there does not exist a fully 
polynomial approximation scheme for the partition function of the hardcore model with fugacity A 
for graphs of maximum degree at most d. 

Phase transitions of spin systems have been known to exactly determine the region of rapid 
mixing in a number of systems including the ferromagnetic Ising model on Z^ [21j and on the d- 
regular tree [3]. The first such example on completely general bounded degree graphs was recently 
established by Mossel and the present author |27j showing rapid mixing of the Glauber dynamics 
of the ferromagnetic Ising model on graphs of maximum degree d when {d — l)tanh/3 < 1. The 
threshold {d — l)tanh/3 = 1 is a statistical physics phase transition, the uniqueness threshold for 
the Ising model on the d-regular tree. 

Slow mixing of MCMC algorithms do not by themselves imply hardness of approximating the 
partition function. Indeed, in the ferromagnetic Ising model the mixing time of local reversible 
Markov chains may be exponential but nonetheless there is an FPRAS by the famous algorithm 
of Jerrum and Sinclair [16j . However, unlike the hardcore model or indeed the anti- ferromagnetic 
Ising model which do exhibit phase transitions, the ground states of the ferromagnetic Ising model 
are trivially found. 

While phase transitions exists on many infinite graphs, it is the uniqueness threshold on the 
tree that appears to determine the onset of computational hardness in general graphs in a number 
of models as they represent the extreme case for correlation decay in graphs for many models. 
Sokal |30j conjectured that uniqueness on the d-regular tree for the hardcore model implies unique- 
ness on any graph of maximum degree d. This conjecture was established in [M] which further 
showed that for any 2-spin system strong spatial mixing on the d-regular tree implies strong spatial 
mixing on all graphs of maximum degree d. Indeed for most, although not all, spin systems the 
regular tree is expected to be the limiting case for extreme correlations amongst all graphs of max- 
imum degree d (see e.g. [31] for more details). The emergence of long range correlations appears 
to be a necessary prerequisite for hardness of sampling and this motivates the conjectures that the 
uniqueness threshold on the tree determines the onset of computational hardness. 

In this paper we establish a form of the second heuristic prediction on a modified random bipartite 
graph. We show that on a polynomial sized set of vertices the spins are close to a product measure, 
conditional on the phase in the L°° distance on measures. Being able to treat large numbers of 
vertices as conditionally independent given the phase plays a key role in our reduction. While some 
results of this nature have been established previously (see e.g. |8l[25|) this is the first example we 
are aware of where the number of conditionally independent sites grows polynomially in the size of 
the graph. 
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1.2. Proof Techniques. Following the approach of [9] and as suggested in [28| we utilize random 
bi-partite graphs as basic gadgets in a hardness reduction. In those papers the basic unit of the 
construction is the random d-regular bipartite graph. To obtain a sharp result we cannot afford to 
add edges to such graphs (creating degree d+1 vertices) so our basic gadgets are bi-partite random 
graphs, most of whose vertices are degree d but with a small number of degree d — 1 vertices which 
are used to connect to other gadgets. 

We begin by constructing a graph G which is a random bipartite graph with n vertices of degree d 
and m' ~ fi^+'^ vertices of degree d—1 where 9, i/j are small positive parameters. We label the sides 
as "plus" and "minus" and edges are chosen according to random matchings of the vertices on the 
two sides. We denote the phase of the configuration (the random independent set) to be plus or 
minus according to the side which has more elements of the set amongst the degree d vertices. 

With U denoting the set of vertices of degree (d—l) we consider the random partition functions 
Z (rj) giving the sum over A''^' over all configurations with phase it and with au = rj where rj G 
{0, 1} . We show that in expectation the KZ (rj) are essentially proportional to the probabilities 
of a product measure on U whose marignals are given by the marginals of extremal Gibbs measures 
for the hardcore model on the {d — l)-ary tree. Our proof requires that this holds approximately 
for the Z (rj) themselves and adopt the second moment approach of [28] including their use of the 
small graph conditioning method [35j. While still involved, by estimating ratios of quantities in our 
model to quantities calculated in [28] we greatly simplify these computations. We are, however, 
still left with the same technical condition as [28] which we describe in the next subsection. 

Even this approximate conditional independence is not sufficient for our reduction. To this end 
we construct a new random graph G by appending (d — l)-ary trees of height ■(/; log(^_]^ n onto U and 
denote the set oi m k, vP roots of the trees as V which are of degree d — 1. Our proof proceeds to 
show that, conditional on the phase, ay is very close to a product measure. We note that appending 
the trees reweights the probabilities on configurations ajj but it does so in a quantifiable way. 

By construction the spins cry are conditionally independent given ajj- Moreover, the statistical 
physics heuristics imply that the configuration of the neighbourhood around ay should be given by 
an extremal semi-translation invariant Gibbs measure on the tree with strong decay of correlation 
from the root to the leaves of the tree. Based on this intuition, we show that after conditioning on 
the phase the probability that ajj has a non-negligible infiuence on ay is doubly exponential small 
in the height. Through this we can establish its distribution with bounds in the L°° norm. This is 
done by bounding the probability that the spins in a distant level infiuence the root using methods 
from the "reconstruction problem on the tree" (see e.g. [2HI32]). 

The random graph G constitutes our gadget. Given a graph H on up to nP''^ vertices we construct 
H*^ by taking a copy of G for each vertex of H. Then for every edge in H we connect n^^'^ vertices 
between each side of V in the corresponding copies of G maintaining the maximum degree d. 
Since the spins in V are almost conditionally independent given the phase we can estimate the 
effect of adding these edges. An easy calculation shows that the most efficient arrangement is to 
have connected gadgets have opposite phases. The hardcore model on H puts most of its mass 
on configurations whose phases are solutions to MAX-CUT on H. Hence, by the equivalence of 
approximate counting and approximate sampling, this gives a randomized reduction to MAX-CUT. 

1.3. Preliminaries. For a finite graph G with edge set E{G) the independent sets are subsets of 
the vertices containing no adjacent vertices or equivalently elements of the set of configurations 

1(G) = {a£ [0, 1]^ : V(u, v) G ^(G), a,a, = 0}. 
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The Hardcore Model is a probability distribution over independent sets of a graph G defined by 



^g{c7) = ^AtT^^"^° "" l-e/(G) (1-1) 



ZgW 

where Zg{X) = X^o-e/fG) "^ ~ X^o-e/fG) A^^'so'^" is a normalizing constant known as the partition 
function and is a weighted counting of the independent sets. When A = 1 the hardcore model is the 
uniform measure on independent sets and Zg{1) is the number of independent sets of the graph. 

The definition of the hardcore model can be extended to infinite graphs by way of the DLR 
condition which essentially says that for every finite set A the configuration on A is given by the 
Gibbs distribution given by a random boundary generated by the measure outside of A. Such a 
measure is called a Gibbs measure and there may be more one or infinitely many such measures 
(see e.g. |12| for more details). When there is exactly one Gibbs measure we say the model has 
uniqueness. Our main result relates the uniqueness threshold on Td, the infinite d-regular tree, to 
the hardness of approximating the partition function on graphs of maximum degree d. 

The hardcore model on T^ undergoes a phase transition at Xc{d) = ^ u_2)d with uniqueness 
when A < Ac and non- uniqueness when A > Ac [Ej. The following picture is described in [28]. 
For every A there exists a unique translation invariant Gibbs measure fi = fi(i,x known as the free 
measure with occupation density p* = ^{(Jp) for p the root of the tree. When A > Ac there also 
exist two semi-translation invariant (that is invariant under parity preserving automorphisms of 
Td) measures ^+ and /i_ whose occupation densities we denote by p'^ = fi+{ap),p~ = fj,-{ap). 
These measures are obtained by conditioning on level 2i (resp. 2^ + 1) of the tree to be completely 
occupied and taking the weak limit as £ — )• cxd. 

It will also be of use to discuss related measures on the infinite {d — l)-ary tree T rooted at p. 
We define analogously the measures /i+ and /i_ obtained by conditioning on level 2£ (resp. 2£ + 1) 
of T(i to be completely occupied and taking the weak limit as £ — )■ oo. We set q^ and q~ to be the 
respective occupation densities q^ = jl-^-{ap), q" = fi^{ap) of the root p. 

The measure p,± and /1± are naturally related as follows. Let t; be a child of p and denote T^ to 
be the subtree of T"' rooted at v. There is a natural identification of T*^ \ T^ with the {d — l)-ary 
tree T'^ and under this identification the measures satisfy 

/i±(cr G •) =M±('7Td\T„ 6 -kf = 0). (1.2) 

In particular since Up = 1 implies o"^ = for an independent set in T*^ it follows that 

Furthermore, standard tree recursions for Gibbs measures (see e.g. [28]) establish that 



,± 



A(l - q'^y-^ 
■" ~ l + A(l-gT)^-i 

and consequently by equation (jl.Sp . 



' A(l-,^)-^ = Af^^^V"- (1-4) 



It is shown in [28l Section 4] and [9l Claim 2.2] that the following hold for A > Ac: 
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^1) The solutions to h{a) = /3, h{(3) = a with (q, [3) £ T = {(a, /3) : a, /3 > 0, a + /3 < 1} where 

l/d 



h{x) 



[l-x) 



1 



X 



X{l-x] 



are exactly (p^,p~), (p , P^) and {p* ,p*)- These densities satisfy p~ < p* < p^ and when 
A I Ac we have that p*,p^,p~^ — )• l/d. 
(2) The points (p'^jP') and {p~ ,p^) are the maxima of <^i(a,/3) in T where 

$1 (a, /3) = (a + /3) log A - a log a - /3 log /3 - (i(l - a - /3) log(l - a - 13) 
+ {d-l) ((1 - a) log(l - a) + (1 - /3) log(l - /3)) . 



1.3.1. Technical Conditions. We now describe the technical condition necessary for our result. The 
function in question is 

f{a,/3,j,6,e)=2{a + /3)logX + H{a) + Hi{j,a) + Hi{a-j,l-a) + H{/3) + Hi{6,P) 



+ Hi{p-5,l-p) + d 



Hiij, 1 - 2/3 + (5) - H{j) + Hiie,l-2f3 + 5-j) 



+ Hi{a-j - e,l3 - 6) - Hi{a-'y,l -'j) + Hi{a-'y,l - 13 -J - e) - Hi{a-j,l - a) 

(1.5) 

where Hi{x, y) = — x(logx — logy) + (x — y)(log(y — x) — log(y)) and H{x) = H{x, 1) and where / 
is defined in the range (a, /3) G T and 



a-7-e>0,/3-(5>0, 1-2/3 + (5-7-e>0. 
which emerges naturally when calculating the second moment of the partition function. 



:i.6) 



Condition 1.2. The technical condition is that there exists a constant X > such that when when 
Ip"*" — /3|, Ip" — a\ < X the function Qa^pil^^^^) = /(ct)/5)7;<Ji^) attains its unique maximum in the 
set ([L6D at the point (7*,(5*,e*) = {a^,l3'^,a{l -a- (3)). 

The following result of [28] establishes Condition 11.21 when Ac < A < \c{d) + e{d). 

Lemma 1.3 ( ^28^ Lemma 6.10, Lemma 5.1]). For each d > 3 there exists x > such that when 
\a — ^1,1/3 — 3I < X then ga,i3{'y,S,£) has a unique maximum at (7*,(5*,e*) where 7* = a'^,6* = 
/32,e* = a(l-a-/3). 

In Section [5] we give a computer assisted proof which establishes Condition 11.21 in the special 
case of A = 1 and d = 6. Two other technical conditions we make use of in the proof are that 



q~^q (d — 1) < 1 and q^ < 



:i.7) 



Both conditions holds in the regions of interest as we have that q^ ,q^ — s- ^jrj when A J, Ac and 
q^ ~ 0.423, q~ ~ 0.056 when A = 1 and d = 6. The first can be shown to hold for all A > Ac with 
a somewhat involved proof while the latter is unnecessary but somewhat simplifies the proof. 
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1.4. Comments and Open Problems. The main open problem, of course, is to remove the 
A < Xc{d)+e{d) condition, ideally with a proof avoiding the second moment analysis. Alternatively, 
one could try and establish the Condition 11.21 for all A > Ac and d> 3. 

Another natural problem is to establish the correspondence between computational hardness and 
phase transitions in the anti-ferromagnetic Ising model. While calculations of the style of [28| are 
not available and are likely to be even more challenging, it may be possible to avoid them. Indeed 
results of [25j already imply conditional local weak convergence of the configuration but not in a 
strong enough form to complete necessary reduction. 

In Section [2. II we detail the construction for G and show how it. In Section[3]we analyze the first 
and second moments of the partition functions Z^ {t]) . In Section U] we analyse the reconstruction 
problem on the tree and establish the conditional distributions of av- Finally in Section [5] we sketch 
the computer assisted proof that Condition 11.21 holds when d = 6 and A = 1. 

1.5. Acknowledgements. A.S. would like to thank Elchanan Mossel for his generous encourage- 
ment, guidance, support and advice with this project and also Dror Weitz for helpful discussions. 
The worked was initiated when the A.S. was a student at UC Berkeley where he was supported by 
NSF CAREER grant DMS-0548249 and by DOD ONR grant (N0014-07-1-05-06) 1300/08. 

2. Proof of Theorem [T] and [2] 

In this section we first describe the construction of our base random graph G which will be the 
basic gadget in our reduction. We state a theorem describing the properties of the hardcore model 
on G and then proceed to show how this establishes the reduction for Theorems [T] and [2l 

2.1. Construction of G. We begin by constructing a random bi-partite (multi)graph G = G{n, 9, Tp) 
where n is a positive integer and < 0, ^ < g are positive constants which will be chosen to depend 
on A and d. This graph will be the basis of our construction of G. 

• The bipartite graph is constructed in two halves which we will call respectively the plus 
half and the minus half each with n + m' vertices where m! = {d— 1)^ °Sd-i"J+2L^ ogd-iJ"^ 

• The vertices of each side are split into two sets W^ and U^ of size n and m' respectively. 
We label the vertices of U^ hy u^ , . . . , u^, . 

• We connect d—1 edges to each vetex by taking d—1 random perfect matchings of W^UU^ 
with W^ U U^ and adding an edge between each pair of matched vertices. 

• We take one more perfect matching of W~^ with W~ and add an edge between each pair of 
matched vertices. 

In this construction the vertices in W = W^ U W~ are of degree d and the vertices in [/ = 
[/"*" U 'W~ are of degree d—1. Note that in this construction there will be multiple edges between 
vertices with asymptotically constant probability bounded away from 1. However, in the hardcore 
model multiple edges are irrelevant and we simply treat them as single edges (some degrees will be 
decreased but this will not affect our proof). 

We now complete our construction of G = G{n, 9, ip) by adjoining trees onto U^ and to U^ . 

• Construct a collection of m = {d — l)L^'°Sd-i"J disconnected {d — l)-ary trees of depth 
2[^ log(^_i n\ rooted at vf,..., v^. The total number of leaves of the trees is m'. 

• Adjoin this collection of trees to U~^ by identifying each vertex of U^ with the leaf of one 
of the trees. Denote the set of roots as V~^ which are vertices of degree d — 1. 
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• Perform the analogous construction on U~ to complete G. 

This construction yields a bi-partite graph of maximum degree d with m vertices of degree d — 1 on 
each side. We now consider a the Hardcore model Pg{'^) on G. Our construction is a modification 
of the model considered in [i28] where they showed that on a. a. a random bi-partite d-regular graphs 
the probability of "balanced" sets is exponentially small. This is also the case for our construction 
and we define the phase of the configuration as 




Y = Y(a) :-- . 

We define the product measure Qy (respectively Q^) on configurations on F = V^ U V^ so that 
the spins are iid Bernoulli with probability q~^ (resp. q~) on V^ and q~ (resp. q'^) on V~ , i.e., 

Q±((jy) := {q^)^vev+ ^- (1 _ q±)^-T.^ev+ ^- (gT)E„6v- ''^ (1 - gT)m-E„ev- ^- . 

We define Qu on U = U^ U U^ similarly. With these definitions we establish the following result 
about hardcore model on G. 

Theorem 2.1. For every d > 3 when Xc{d) < A and when Condition \1.2\ and equation ()1.7p hold 
there exists constants 9{X,d),ip{X,d) > such that the graph G{n,9,ip) has (2-|-o(l))n vertices and 
satisfies the following with high probability: 

• The phases occur with roughly balanced probability so that 

Pg{Y = +)> -,Fg{Y = -)>-. (2.1) 

n n 

• The conditional distribution of the configuration on V satisfies 



max 

(TV 



¥G{av\Y = ±) _^ 



Qvi^'v) 



< n-"^^. (2.2) 



The proof of this theorem is deferred to Section [H 

2.2. Reduction to Max-Cut. We now demonstrate how to use Theorem 12.11 to establish a re- 
duction from sampling from the hardcore model to Max-Cut. Let if be a graph on up to -^^n^'^ 
vertices. With a random bi-partite graph G = G{n, 6, ip) constructed as above we define H^ as 
follows. 

• Take the graph comprising \H\ disconnected copies of G and identify each copy with with 
a vertex in H labeling the copies {Gx)x&H- Denote this graph by H^ . We let V^ and V~ 
denote the vertices of Gx corresponding to V^ and V~ . 

• For every edge (x, y) in the graph H add n^^' ^ edges between V^ and V^ and similarly 
add n^^'^ edges between and V' and V~ . This can be done deterministically in such a way 
that no vertex in H has its degree increased by more than 1. Denote the resulting graph 
by iJ'^. 

The resulting graph has maximum degree d. For each x £ H we let Yx = Yx{a) denote the phase 
of a configuration a on Gx- Let y = {Yx)xeH £ {0, 1}^ denote the vector of phases of the Gx- 
Denote the partition function given the phase y by 

ZHo{y')= Yl AHi(3;(a)=y). 
(Tei(HG) 
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Lemma 2.2. Suppose that G satisfies equations ()2.ip and ()2.2p of Theorem \2.1\ Then 



Zsaiy') 



Z 



Pg(y = +)S^eff ^y^=+ . Pg(^ = -)^^sH ly^=- > 



n 



„e/4 



hg 



and 



ZHo{y) 



{Ch + o{1)) 

,36/4, 



:l-g^ 



(l_(^+)2)(l_(^-)2) 



n^»/4Cut(y) 



(2.3) 



(2.4) 



znoiy) 

where Ch = [(1 - (g+)')(l - (g-)^)]""''''^''^ and where Cut(y) = #{(x,y) G E{H) : X ^ J^^} 
denotes the number of edges in cut of H induced by y' . 

Proof. Since the graph H'^ consists of a cohection of disconnected copies of G, the distribution 
of a configuration on H is given by the product measure of configurations on the {Gx)xeH- In 
particular the phases are independent and so 

zsoiy') 



'HG 



V {y{(T) = y) = ^g{Y = +)^-sH ly^=+ . p^(y = _-)^.eH ly^=- > 



n 



,e/4 



which estabhshes equation (j2.3p . Now the ratio of the partition functions in (|2.4p is exactly the 
probability that the configuration a sampled under IP^g is also an independent set for H'-' after 
adding in the extra edges, that is 

|^ = p^.(.G/(//^)lw = y) 

= Ffja (y{v,v') G EiH^) \ E{H^),a,(j„, / 1 | y{a) = y 



Now by equation (j2.2p . conditional on the phase y' the spins of cru^g^vi ai"e asymptotically condi- 
tionally independent with probabilities q^ or q~ depending on the phase. It follows that 



^HG [^{v^v') G E(F^) \ E{H''),a,a,> / 1 | y{a) = y' 
= (l + o(l)) J] f'HG(.cT,c7,,^l\y{a)=y) 

If {x,x') G E{H) then by direction calculations and equation (j2.2p 

n PjjG(^.^.'/iiK^) = y) 

i;eG:,,-u'eG^/:{v,-u')e^{^'^) 



(l+0(n-^))((l-(g+)2)(l-(g-)2)) 
(l + 0(n-^))((l-,+g-)2)""'' 

Combining the above estimates we have that 



,39/4 



ZHo{y) 
z^a{y) 

which completes the proof. 



(Ch + o(1)) 



\iYx = y.'. 

iiYx^Yx>. 



{I - {q+?){l - {q-Y) 



Given the previous lemma we now show how to produce the randomized reduction to Max-Cut 
establishing Theorems [1] and [2j 



COMPUTATIONAL TRANSITION AT THE UNIQUENESS THRESHOLD 



11 



TheoremlJl andl^ Let H he a graph on at most ^^n^ vertices. Take an instance of a random 
graph G = G{n, 9, ijj) according to the construction in SectionEUJ By Theorem 12.11 with probabihty 
tending to 1 the graph satisfies equations (|2.ip and (j2.2p . Assume that it does and construct the 



graph H which has at most 0{n^^ ) vertices and maximum degree d. 

Now suppose there exists an FPRAS for the partition function for the hardcore model with 
fugacity A on graphs of maximum degree d. We now use the equivalence of approximating the 
partition function and approximately sampling for the hardcore model described in the introduction. 
In polynomial time we may approximately sample from the hardcore model on H'^ to within 
6 of the Gibbs distribution in total-variation distance for any 6 > 0. Let a' denote such an 
approximate sample. We may couple a' and with a distributed according to the Gibbs measure so 
that F{a' j^a) <d. We now consider the phase of a. Let y,y" £ {0, 1}^ such that 

cut(y) > cut(y'). 



Then by Lemma 12.21 we have that 



r{y{a)=y") Zfjciy) - 



ZHciy") 

>(l + o(l))n-" 



e/4 



{l-q+q-f 



(1 - q+q- 



^)(1 



-,n^o/-i[Cut{y')^Cut{y")] 



(l-(g+)2)(l-(g-)2) 



n3»/4[Cut(y)-Cut(y")] 



(2.5) 



As we have that < g" < q+ < 1 if follows that (l-g+g-)2-(l-(g+)2)(l-(q-)2) 
and hence 



{q+-q-y>0 



(1 - q+q- 



(1 - (g+)2)(l - (g-)2) 



> L 



Therefore, for large enough n by equation (j2.5p it follows that 



nyja) = y) 

F{y{a) = y") 



> 



(1 



q q 



-^2 



(l-(g+)2)(l-(g-)2) 



,39/4 



>4^ 



,e/4 



Since the size of {0, 1}I^I is only 2""'''' it follows that with probability at least l-2l^l that Cut(3^(cr)) 
attains the maximum value. Hence with probability at least 1 — 5 — o(l) the phases y{cr') of the 
approximate sample a' also attains a maximum cut in H. As such we have constructed a randomized 
polynomial-time reduction from approximating the partition function of the hardcore model to 
constructing a maximum cut. It follows that unless RP=NP there is no polynomial-time algorithm 
for approximating the partition function of the hardcore model for Xc{d) < A < Xc{d) + e{d) on 
graphs of maximum degree d or when A = 1 on graphs of maximum degree 6 or more. ■ 
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3. The partition function of G 

In this section we analyse the hardcore model on the random bi-partite graph G and in particular 
consider the effect of conditioning on the spins in U = U^ U U~ . For r/ E {0, 1} we define ZQ{r]) 
to be the partition function over configurations whose restriction to C/ is r/, that is 

Our analysis borrows heavily on hard computations carried out in |28] . There they considered a 
random d-regular bipartite graph where each side has n vertices and the edges are chosen according 
to d independent perfect matchings of the vertices of the sides. They denote Z°-'P to be the weighted 
sum over configurations of the graph with an and fin vertices on the plus and minus sides of the 
configuration (for a,/3 such that an, fin are integers). We will denote their quantity by Z^'^-^^. In 
the same spirit define 

Lemma 3.1. For any (a,/3) in the interior of T and all t] G {0, 1}^ we have that: 



EZj''(r/) = (l+0(n 



-1/2 



))C* A 



l-a-/3 
1-/3 



A 



l-a-/3 



1 



a 



EZ, 



a,(3 
MWW 



(3.1) 



where 



C* 



(l-«)(l-/3) 
/3 



1 



a 



and where Vj denotes X^„g(7± ?/«• 



Proof. We follow the approach of [28] in estimating these quantities. In total there are (^) ( IJ^) 
choices of configurations on W with an sites on the top and fin sites on the bottom. Then by 
calculating the probability that a perfect matching does not connect two I's of the configuration 
we have that 



'f' 



EZT{ri) 



\an+l3n+rj'^ +ti 



n 



an 



n 

fin 



/n+m'—l3n—ri 
\ an+ri+ 



/ n+m' 



while by [28] we have that 



d-i 



tn-j3n\ 
\ an ) 

(") 
\anJ 



(3.2) 



^MWW 



\an+(in( ^ 



an 



n 
/3n 



\ an ) 

(") 
\anJ 



Now since \U\ = 0{n^''^) it follows from Lemma 13.21 below that 

d-i\ T 






(l + 0(n-i/2))(^* IX 



l-a- f5 
1-/3 



l-a- fi 



a 



d-i\ ''^ 



To complete Lemma 13.11 we give the following lemma which a simple expansion of factorials 
which will use repeatedly throughout this section. 
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Lemma 3.2. When < b < a are integers and x^ + y^ < min{6, a — b} then 



/a+x\ , , 9 9\\/ \ X / I 

V6+J/J (-, , r^l x^'+y' W a \ fa-b 



Kb+y) 



,, y yniin{6, a — 6} yy \a — b 

Proof. By expanding out factorials we have that 

iZy) _ ia + x)\ b\ (a -by. 

Q a! {b + yy.{a-b + x-y)l 



^ + --ttzHj Z—. M- ■ (3-3) 



y 



x—y 



n(i+i)j^r.n(i-j)j^(«-^-n(i+^)^ 

X + y W ( fl \ I a — b 



1 + 



niin{6, a — b} J J \a — b 



We now sum over (a, /?) and define the conditional partition functions as 

andZ|=E,^|(r?). 

Lemma 3.3. For every d > 3 there exists constants 0*{X, d),'ip*{X, d) > such that when Xc{d) < A 
and < 6{X,d) < 6*{X,d),0 < ip{X,d) < 'ijj*{X,d) then the expected partition functions satisfy 



sup 






oil) (3.4) 



and 



EZ+ = (l + o(l))EZT. (3.5) 



Proof. Recall from Lemma 13.11 that 

EZ-/ir,) = il + Oin-y^))C*^xl^^^^y~y (^A (i^^)"'^ EZ-^^ (3.6) 

and that by [281 Proposition 3.1], 

^^MWW ^ exp($i(a,/3)n) 

where the approximation holds up to a polynomial factor in n. In the proof of [9l Claim 2.2] it 
is shown that for fixed a (resp. /3) ^i is maximized by setting /? = h{a) (resp. a = h{f3)) where 
h{x) = (1 - x)[l - (x/(A(l - x)))^/"^] was defined in the Section[I31 Recall that in {(a,/3) G T : 
o: > f3} the function <l>i is maximized at {p'^,p~). Clearly we have that the functions $i(a,/i(a)) 
(resp. $i(/i(/3),/3)) are analytic in a (resp. /3) when in a neighbourhood of p"*" (resp. p~). It follows 
by expanding as a Taylor series and noting that {p'^,p~) is a local maxima that for some integer 
i >2 and constants C, e > we have that 

\<^i{a,hia))-<^i{p+,p')\ >C\a-p+\' 
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when \a — p'^\ < e and similarly for <I>i(/i(/3),/3). This of course implies that when ||(q,/3) — 
(p+,_p")||oo < e then 

|c&i(a,/?) - $i(p+,p-)| > C||(a,/3) - (p+,p-)||L- 

Hence it follows that for large n, 

^ Zf{r^) < exp (-f n^/^) EZt(r?). (3.7) 

a>/3,l!(a,/3)-(p+,p-)||oo>n"5I 

Setting 0*(A, d) = '4'*{X, d) = -^ we have that \U\ < nsi and hence 



1 — ay / \\1— p+ 

Combining equations (|3.6p . (|3.7p and (|3.8p we establish that 

Ezg(,) = (1 + o(i))c* [a ( ^ 7_";/' )"T (a( '7_";/' )"T ez^^w 



= (l + o(l))C7*^^^j Vr^j ^^^^^ 

~ (1 _ g+)m'(l _ q-)m' ^t/W^^MWW 1-^-^^ 

where EZj^^^ denotes J2a>i3^^MWW^ where the second line follows by equation p.3p and the 
final equality holds by the definition of Q'^. Hence we have that 

EZ+ - ^^ + "^^^^^* EZ+ 

^^G ~ (1 - g+)™'(l - g-)m'^^MWW 



and that 



sup 



^'>' QS(,) 



EZ± ^^' 

(sr 



oil) 



The analogous bound holds for Z- (rj) which establishes equation ([3 

Note that because of the slight asymmetry in the definition of Z'^ and Z'- equation (j3.5p is not 
immediate by symmetry. It follows from the fact that when A > Ac by symmetry we have that 

y"EZ"'"=EZt-EZT 

a 

and 



7 + 



J^EZ^'" < exp(-0(n))EZi 
since the maxima of <I>i(a, /3) is not achieved with a = /3. This completes the lemma. 
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3.1. Second Moment Analysis. We now proceed to analyze the second moment of the partition 
function. In [28j they showed that the second moment is given by 



E 



^MWW 



iv) 



n 



;^2(a+,3) f 



n 

/3n 



7,5 



an 



"fn 



(1 



a 



a)n 



(3n 
6n 



(1 - f3)n 
il3 - 6)n 



fil-2l3+S)n 



\7n 



) 



'e' 



(l-2/3+5-7)nW {I3-S)n \ 
en / V{o— 7— £)n/ 



|'(l-/3-7-e)n\ 
V (o— 7)n / 



|'(l-7)n\ 
v(a— 7)n/ 



^(l-a)n\ 
\(a—'y)nl 



(3.10) 



where the sums run over 7, 5, e such that 771, 5n, en and equation (jl.6|) holds. Equation (j3.10p 
should be interpreted as follows: The first line represents the number of ways of choosing a pair 
of configurations, both with size a on the plus side and /3 on the minus side with overlaps of 7 
on the plus side and 6 on the minus side. The second line gives the probability that the pair of 
configurations are both independent sets in the random graph (see [28J for the interpretation of the 
sum) . Here the role of Condition 11.21 comes into play. A simple approximation gives that 



exp 



(n/(a,/3,7,5,e))«A2("+/^)f" 

\an 



n 

(3n 



an 



7n 



(1 - a)n\ //3n\ /(I - /3)n 
(a - j)n) \5n) V(/3 - 6)n 

d 



^(l-2/3+^)n^^ ({l-2l3+5-y)n\ f W-S)n \ /(l-/3-7-e)n 
\ en J \{a—"/—e)n/ \ (0—7)71 



■yn 



rn^ 
y^n/ 



/(l-7)n\ 
V(a— 7)n/ 



i{l-a)n\ 
\{a—'y)n) 



up to polynomial terms in n. As such the maximum of / plays a crucial role in the second moment 
analysis. The following result is by [28^ Lemma 3.3]. While they only stated their result for (a, /3) 
close to (1/d, 1/d) it is easy to verify that their proof holds in a neighbourhood of {p~ ^p'^) whenever 
Condition 11.21 holds. 

Lemma 3.4 ([28;, Lemma 3.3]). For each d > 3 suppose that Condition \1.S\ holds. Then there 
exists some x > such that when |a — p~|, |/3 — p^| < x we have that, 



E z; 



o,/3 
MWW 



'EZ 



a,/3 
MWW 



where 



.a,/3 



^a-.P 



(1 - a - /3 - afiy 



[{I- a- 13 + 2a/3)(l - a - /3)] V [(1 - a - /3 + da/3)(l - a - /3 - {d - 2)a/3)] 2 



Next we show the analogous result for G conditioned on au by estimating the ratio of the second 
moments of the partition functions of the graphs. 

Lemma 3.5. For each d>3 suppose that Condition \1.2\ holds. Then there exists some x > such 
that when \a — p~\,\P — p'^\ < x we have that for all rj G {0, 1} , 



E(Z|'^(r?) 



E z; 



^MWW 



(l + o(l))(C*)MA 



l-a- P 
1-/3 



d-V 



27?- 



l-a- P 
1 — a 



d-l 



2r]-^ 



(3.11) 
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and hence 



Proof. Repeating the analysis of [28] the analogous 

7,5 



T°'^. 



an\ /(I - a)n\ //3n\ /(I - f3)n 
jnj \{a — 'j)nj \6nJ \{j3 — 6)n 



({l-2l3+S)n+m-r]^\ /(l-2l3+S-^)n+m-ri+ -r]'\ / (P-5)n \ /(l-,3-7-£)n+m-r;+\ "" 



7n+»?+ I V~^ V en / \(a— 7— £)n/ V (a— 7)n / 

/ n+m \ Z_^ /(i_^)„_|_m-r;+-r?-\ l{\~a)n+m-r]-\ 

W"+f?+^ e I, (a-7)n / V (a-7)n / 



(3.12) 



By Lemma [TTBl the unique maxima of ga,pk^^^^^) is (t*,'^*;^*) and it was shown in [25] that ^q,^^ 
decays quadratically from this point. Consequently, with 

A = {(7, 6,e):h-j*l\S-6*\,\e-e*\< n-'/^} 
the contribution from terms with (7, 5, e) ^ A is exp{—Q{n^''^)) and so can be omitted. Setting 

/{l—2(5+S)n+m—ri~\ /(I— 2/3+5— 7)n+m-»y+— »;~\ / {fS-S)n \ /(l-/3-7-£)n+m-r?+\ 



^"'/^(ry) = A2("+^)- 



7n+r7+ / V en / \{o— 7— e)n/ V (a— 7)™ 

/ n+m N /{l-'y)n+m-ri+ -r]~ \ /{l—a)n+m—-q-\ 

{■yn+r]+) [ (a-'y)n ) \ (a-'y)n ) 



K 



{a—'y)n J \ {a—'y)n 

n\ /{l-2(3+5-y)n\/ (I3-S)n \ /(l-/3-7-£ 
a, 13 I 7n / V en ) \{a—'y—s)n/ \ {a—j)n 



MWW /n\ /(i_^)„N /{l-a)n\ 

W"'' l(a-7)n/ \{a-y)n) 

and recalling that 7* = a^, 6* = /3^, e* = a(l — q — /3) we have that for (7, 6,s) ^ A 



^6^^"^^ _n , .n^^^2(v++v-) f 1-2/3 + 5 



(1 + 0(1)) A2(''+ +''')/' 



m—rj 



k"'^ \l -213 + 6--/ 

l_2/3 + 5-7V^ / ^ V+^ l_2/3 + 5-7 ^'"-''^''" 

(1-7) 



7 / vl-7/ \l-2l3 + 6 -J- e 

1 — p— 7 — e\ ' ' f I — a\ ' /I — 2a + 7^ ' 



1 — /3 — a — ey \1 — 7/ \ 1 — a 



where the first line follows by equation (j3.3p and the second follows by approximating (7, 5, e) with 
(7* , 5* , e* ) and simplifying. Now comparing equations (j3.10p and (j3.12p (noting here that we can 
neglect terms with (7, b^ e) A) we have that 

which establishes equation (j3.1ip . Equation ()3.12p then follows from Lemmas 13.11 and [37 
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3.2. Small Graph Conditioning Method. As we showed in the previous section the ratio of 
second moment and the first moment squared of Z'^ (rj) converges to a constant r"'^. Unfortu- 
nately, since r"'^ > 1 we can not directly apply the second moment method to get high probability 
bounds. Instead we follow the approach of |28j and use the small graph conditioning method. Our 
proofs differ very minimally from theirs and as such we only comment on the necessary changes. At 
a high level, the method says that small cycles in the graph "explain" the variance which provides 
good lower bounds on Z'^ [r]). The following is taken from [28j which itself is presented as a special 
case of a results of [3S] and [K] . 

Theorem 3.6 ([^S", Theorem 7.1]). Let Aj > and 6i > be real numbers for i = 1,2, .. . Let 
uj{n) — )• and suppose that for each n there are random variables Xi = Xi{n),i = 1,2, .. . and 
Y = Y{n), all defined on the same probability space G = Gn such that Xi are nonnegative integer 
valued, Y is nonnegative and EY > (for n sufficiently large). Suppose furthermore that 

(1) For each k > 1, the variables Xi, . . . ,X^. are asymptotically independent Poisson random 
variables with EXj — )• Aj 

(2) For every finite sequence mi , ■ ■ ■ , m^ of nonnegative integers, 



EY 



n(A.(l + <5.)r (3.14) 



where [X]m = YYlLo (^ ~ ^)' denotes the falling factorial. 

(3) That Y, Xi6f < oo, 

(4) That EY^/{EYf < exp(^ XiSf) + o(l) as n ^ oo. 

Then ¥{Y > uj{n)EY) -^ 1. 

We set Y = A^*-"^'^^"Z'? (r/) and let Xi be the number of cycles of length i whose vertices lie in 
W (which is of course when i is odd) . The following lemma has an essentially identical proof to 
the proof of [281 Lemma 7.3] and follows from standard methods. 

Lemma 3.7. For even i the number of cycles are asymptotically Poisson with means Aj = r{d,i)/i 
where r(d,i) counts the number of proper d-colourings of a cycle of size i. 

Next, the main step is to determine the limit of equation (I3.14p . The proof of the following 
lemma follows with only very minor modifications form that of [28, Lemma 7.4 and 7.5]. 

Lemma 3.8. For all (a, /3) in the interior of T and all rj G {0, 1} we have that: 



i=l 



where Xi = r{d,i)/i and 



,(l-a)(l-/3) 
for even i. 

Proof. We follow as far as possible the proof of [281 Lemma 7.4]. We consider just the case where 
a single rm = 1 and the others are all as the extension to general mj's is exactly as in |28i Lemma 
7.5]. We follow the notation is from |28] with only slight modifications to our setting 
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• Let T G {0, 1}^ with Y.v,&w+ '^^ = ""■ ^^^ with Y.w&v- ^^ = f^^- 

• Denote by ^ a proper d-edge-coloured rooted, oriented i-cycle (from amongst the r{d,i) 
possibihties) , in which the vertices are 2-coloured, black and white, with no two black 
vertices adjacent. The color of the edges will prescribe which of the d perfect matchings 
an edge of a (potential) cycle will belong to. The black vertices will prescribe which of the 
cycle vertices are members of {w £ W : T^^, = 1}. 

• We let C denotes a position that an i-cycle can be in (i.e. the exact vertices of W it traverses, 
in order) such that the prescription of the vertex colors of ^ is satisfied. (Note this was 
denoted as tj in [28]). 

• Denote by PI is the probability that a random graph G contains a cycle C in the given 
position C with the edge colors prescribed by ^ in accordance with which matchings contain 
the edges of C. 

• We denote by P2 the conditional probability that in the random graph G that the set 
{w £ W : Tw = 1} U {u G [/ : 7/u = 1} is an independent set, given that it contains C as in 
the definition of PI. 

• Denote by P3 the probability that in the random graph G the set {w £ W : T^ = l}U{u £ 
U : r]u = 1} is an independent set. 

Analogously to equation (18) of |28] we have that 

EYXj _ 1^^ P1P2 
EY i ^^ P3 ^ ' ' 

as the probabilities are independent of T. It is immediate from the definition that PI = (1 + 
o{l))n-\ 

Now closely following the notation of |28j for k = 1, . . . ,d let e{k) denote the number of edges 
of colour k in ^. Denote by f±{k) the number of black vertices adjacent to edges of colour k in 
the sets {w £ W^ : T^ = 1}. Assuming that ^ is compatible with {w £ W^ : T^ = 1} then the 
probability that the remaining edges also respect the independents sets is given by, 



ut\ 



1 (n+m'-l3n-r]- -e{k)+f-{k)\\ fn-l3n-e{d)+f-{d)-^ 
an+ri+-f+{k) )) \ an-f+{d) 



P2= ^ ' — ■' ^^^-' -/ '^^j^y:^ _^ (3_lg) 

/t-to!— 1 / n+m'-e(k) \\ ( n-e(d) \ 
\\.\.k=l \an+ri+-f+{k)) J \an~f+{d)) 

Now by Lemma 13.21 we have that 

I an+ri+-f+ik) ^ _ (^ _^ ^^ )) ' n + m - f3n - 7] \ 



(n+m'-i3n-'n-\ ^ " K^ + m' - Pu - r]~ - an - t] 

n + m — pn — T] — an — 7]^ \ 



an + r/+ 
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where we used the fact that m,r/"*",7/ = 0{n '^'^) and e{k), f^{k), fj^{k) = 0(1). Similarly we 
have that 

-e{d)+f-id) /I _ fl _ ^A -f+{d) 



fn-Pn-e(d)+f.{d)\ 
\ an-f+{d) I 
fn—l3n\ 
\ an I 
I n+m' — e{k) \ 
\an+r]+-f+{k)) 
/ n+m' \ 
\an+r}+/ 

( n-e{d) \ 
\an+-f+{d)) 



\an/ 

By equation (|3.2p we have that 



(l + o(l)) 
(l + o(l)) 
(l + o(l)) 



1-/? 
I- 13-a 

I \ -e(fe) 



I- 13-a 



a 



1 — a 



1 — a 



-e{d) 



1 — a 



a 



1 — a 



a 



-f+{k) 



'f+id) 



P3 



In+m'—fin—rj 
\ an+ri+ 
I n+m' \ 
\an+ri+ ) 



d-1 



/n—l3n\ 
\ an J 

(") 
\an/ 



(3.18) 



(3.19) 



1 v^rf 



Now let j±(0 = 2 Sfc=i /±(^) denote the number of black vertices in V according to ^ and recall 
that i = Y1=i e(^)- Combining equations ([3l^ . (l3T7D . (|3l^ and (|3l^ we have that 

P2^ (l-a-/3)^-2.-2+ 

P2, ^ ^ ^ '' (1 - a)^-2i+(l - /3)i-2i- ' 
Now letting p/;^^ww (^g^ote the corresponding probabilities in Lemma 7.4 of [28j we note that 



Hence 



PI = (1 + o(l))Pl^ww 



P2 
P3 



(l + o(l)) 



P2 



MWW 



P3MWW 



Ey 



(i+^WJiEE 



« C 



PI P2 
P3 

I MWW DoMWW 



= (1 + «(!))- EE pgMWW 

« C 

= {l + o{l))\i{l + 5i), 

where the final term is the main result of [28, Lemma 7.4]. The complete result for general mj's 
follows similarly to [281 Lemma 7.5] which completes the lemma. I 

Lemma 3.9. If d > 3 and (a,/3) is in the interior ofT and the function ga^p achieves its unique 
maxima in (II. Gh at {a^ ,0^ ,a{l — a — /3)) then for all rj G {0, 1} , 



7a,/3/ 



7a,l3/ 



supP(^Z|'^(7?)<— EZ|''(7?) 



0. 



(3.20) 



Proof. The result follows from an application of Theorem l3.6l taking the i to be even, Aj = r(d, i)/i 

Condition (1) of the theorem holds by Lemma 13.71 Condition (2) holds 



and 6j 



a/3 



j/2 



^(l-a)(l-/3)^ 

by Lemma 13.81 Conditions (3) and (4) hold as a consequence of Lemma 13.51 and [28, Lemma 7.6]. 
Taking uj{n) = -y=, the result follows. ■ 

In Lemma 13.31 we gave estimates of the expected conditional partition functions. In this sub- 
section we of the small graph conditioning method results and give with high probability type 
estimates for the conditional partition functions. 
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Theorem 3.10. For every d > 3 and A > Ac such that Condition \l.S\ holds there exists a positive 
constant e{d) > and constants 9*{\,d),'ip*{X,d) > such that the partition functions satisfy the 
following asymptotic almost sure statements, 

sup P f Z±(r?) < -^EZp,iij)) ^ 0. (3.21) 

Proof. Condition 11.21 guarantees that for (q,/3) in a neighborhood of (p~,p^) that ga,i3 achieves its 
unique maxima in (jl.6p at (a^,/3^,a(l — a — /3)). For sufficiently smah 5 > then by Lemma [37 
we have that 

sup sup P (z%'^{t]) < -^EZ%^^{r,)] -^ 



V (a,/3)65 

where S = {(a,/3) : ||(a,/3) — {p'^ ,p~)\\oo < ^} and hence 



^ \(a,/3)ecS ^ (a,/3)& 



supP I > : Z^^^irj) <^ y: EZ^'^ir^) | ^ 0. 



By equation (|3.7p we have that for all rj £ {0,1} 



and hence we have that 



supP TzJCt?) < -^EZ±{ri)] ^ 0. 



The analogous statement for Z- (77) holds similarly which completes the lemma. I 

4. Reconstruction on the tree 

Our proof now takes a detour through the reconstruction problem on the tree. This problem 
concerns determining which Gibbs measures on the tree are extremal, or equivalently when the 
tail cr-algebra is trivial or when point-to-set correlations converge to in the distance of the point 
to the set p6]. In our setting the measures fL± are extremal so we automatically have that non- 
reconstruction holds. We will use facts about the rate of decay of point-to-set correlations to 
establish that ay is essentially independent of au conditioned on the phase. In most cases the 
reconstruction problem has been considered in the case of the translation invariant free measure 
(see [5] for recent progress on the hardcore model) but we will be interested in the case of the semi- 
translation invariant measures fl± on T^ and as such results from the literature do not directly 
apply here. 

The reconstruction problem has for the most part been studied in the case of Markov models on 
trees with a single transition kernel M. In this theory the key role is played by the X^, the second 
eigenvalue of the transition matrix. The famous Kesten-Stigum bound [18y26j states that there is 
reconstruction when A^(ri — 1) > 1 while results of [T^ show that if non-reconstruction holds and 
Xl{d—1) < 1 then point to set correlations decay exponentially quickly. In our setting, however, the 
Gibbs measure is semi-translation invariant and the Markov model is given by a pair of alternating 
Markov transition kernels, M defined below. 

With minor modifications the proof of [7] (or also [32] or |14j ) can be adapted to the semi- 
translation invariant setting. Here the role of A* is played by the second eigenvalue of M1M2 and 
there is reconstruction when Xl{d — 1)^ > 1 and exponential decay of correlations when there is 
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non-reconstruction and A^((i — 1)^ < 1. The term (d — 1)^ is explained by the fact that this this 
the branching from two levels of the tree. Using the methods of [32] which build on the work of [7] 
we establish the necessary decay of correlations result. 

While we will be interested in the measure jl it will be most convenient to work first on an 
adjusted Markov model ^^ on the tree T^ taking values in {0, l}"""'' and then transfer results to jl. 



The spin ^^ is chosen according to 



■± 



1 



1 



4^ = 1 



p 



,± 



For the other vertices of the graph the values of ^^ will be propagated along the tree given though 
Markov transitions given by alternating transition kernels. Specifically if vertex u is the parent of 
V in the tree then the spin at v is defined according to the probabilities 

p{C = 3\'yu = i) = MtiSl 

for s G { — 1,+!} and i,j G 0, 1 and where \v\ = d{p,v) and 



M+i 



1 







M-^ 



1-q- 
1 







Viewing T^ C T^ we have that the measure of ^ is simply the projection of fi to Td (had we 



instead chosen P 



t 



1 



q the ^ would be given by fl ). It follows that 






if s(-l)l''l =+1 
if s(-l)l"l = -1 



For a vertex u G T let T^ denote the subtree of descendants of v (including v). Observe that the 
measure ^'^ restricted to T^ is equal in distribution to ^*(~^) " on T^ appropriately shifted. Now let 
Sy^i denote the set of vertices in T^^ which are £ levels below v and let ^^ ^ denote the configuration 
on Sv^£. For a configuration A on Sy/ and s G {—1,+!} define the posterior function h^ ^ as 

We set 

for s G {+, — }. Now since the measures fi^ are extremal it follows (see e.g. j24j ) that 

X.As'^f^-'^'^'- (4.1) 

Moreover, if ui, . . . , Ud-i are the children of p then by standard tree recursions for Gibbs measures, 



X, 



p,e,s 



p'UUa-p-')-'[^ - Xua-i,s] + (1 -p^) nti f^^«./-i,s + &[i - ^«./-i,s] 



p n»=i[i — 



i-p- 



p'Utill - """-lyl + (1 -P^)nti fl + [^«./-M -p-1(T3^Ffc^ 



A 



(4.2) 



where the second inequality follows from equation (jl.3p . Next, similarly to [32], we set 



Xi^s = IE Xp^i^s -P^- 
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We will let E^ (resp. E*^) denote the expectation conditional on the ^p = 1 (resp. S,p = 0). With 
ui, . . . , Ud-i the children of p we have the following relationships of the X„ ^ ,,. 

Lemma 4.1. The following hold: 

• Conditional on S,p the Xu^/^s o'^e conditionally independent. 






-Xf^ 



Also¥}{Xu,,i,s-p 

We have that xi^g = (p'*)~^]E(Xp^£_s -p'")^- 

For all integers k > 1 we have that 



Ei(X,^,,,,-p-^)'= = 0(x,,_,). 

Proof. The first part follows from the Markov property of ^. The second follows from the fact that 
EX„,_^_i ,, — p^^ = 0. The third part follows from the proof of Lemma 2.2 of [32j. Finally for the 
forth part we have that 



— s\k 



E\Xu,,e^s-p 



Now since \Xpi -s — P *|<lwe have that 



E" 



s\k 



E" 



{Xp,e,^s-p-l'' 



{Xp 



<W 



P 



(X 



-s\2 



p,t,-s 



for k > 3. When A; = 1 we have 



E^ 



(Xp,^ _s - P" 



P -E' 



1 — p 



{Xp,e,-s-P 1 =0{xe^_s) 



while when k = 2 we have that 



E" 



{Xp^e,-s-p-r <{l-p-T'^ {Xp,e,-s-p-r =0(x,,_,) 



which completes the proof. 



'-=A-A{13-l) + {B- if 4 <A-A{B-l) + {B- if 



We now expand out equation (j4.2p as 

A 

B ■" ■"'" ^' ' '- ^' B 
since A < B. We may expand out B — 1 and can express in the form 

d-i 

B-l= Yl Cc^UiXu^-hs - P'T' 

ae{o,i}'*-i «=i 



d-1 



where for some constants Cq. Moreover, c^ = if |a| E {0, 1} where |a| = X]i=i '^i- Since 



(4.3) 



iE'n(^«« 



.p-T'=o(xi"J,_j 



i=l 



it follows from Lemma |4. II that 

E^ [-AiB -1) + {B- if] = 0(xti^ 



(4.4) 
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Similarly we have that 



/-^ 1 — n * 






^p' + 7rr^^''^"''"' + '^^'^'-i'-^^- (4.5) 



Combining equations ()4.3p . ()4.4p and (j4.5p we have that 



5/,s — P ^ 
and after iterating we have that 






,1^ . {d-lf{p^p-^f 



(d - l)\q+q~fx^.2,s + 0(3:i2,s)- (4-6) 



Now by equation ()4.ip we have that xn^s — ?• as ^ — )• and hence by equation (j4.6p it converges 
exponentially fast to as one of our initial assumptions in equation (II. 7p was that q^q~ {d—1) < 1. 
It follows that by the second part of Lemma |4. II that there exist constants Ci{X,d),C2{X,d) > 
such that, 

E|X,,,,, - p^-'A |2 < c^ eM-C2i). (4.7) 

We now define ^^''" for s G {—1, +1} and v G T^ as the Markov model on the subtree T^ with the 
same transition matrices but so that the initial distribution at v is given by 

P [^^.^ = 1] = 1 - P [^^.^ = 1] = q<-^A . 

With this initial distribution £^^'^ is distributed according to the extremal hardcore measure fi^^~^' " 
on T"^. Analogously to ^, for a configuration A on S^^g and s G {—1,+!} we define the posterior 
function h^ ^ as 

for s G {+, — }. By the definition of conditional probability and the Markov property of ^ and ^, 

mr = MC} = A) _ nc} = mr = i) p(eg-- = i) 
nev"" = o\cj = A) p(e;'; = ^i^^'^' = o) " p(e.^'^ = o) 

_ Fjii, = A\if, = 1) rjcr = 1) 

niii = ^le = 0) ■ mv'" = 0) 

^ p(lg = i\ii, = A) rjcr = imif, = q) 

P(e = 0|f;, = A) ■ FiCv" = 0)m^ = 1) ' 
It follows that for some constant C and any configuration A on S^/ that 

KM)-q'^-'^'^'\ < c\K^M)-f^~'A\_ (4.8) 

We define 
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Note that we are taking the posterior function for v but the Markov model rooted at p which is a 
standard object in the recursive analysis of Gibbs measures on trees. By the Markov property and 
equation ()4.8p we have that 



EIX, 



s(-l)l"l|2 



< CK 



\Xr, 



,Ls - P 



;(-l)l"l|2 






s,p 



1 



+ C7E 



<2C(min{/(-i)'"',l-/(-i)""}) E|X 



-1 



s{-l)l"l|2 



i7 = 



^vl,s — P 



and hence we may conclude that there exists constants, Ci,C2 > such that for all v G T^ and 
s G {—1, +1} we have that 



E|X,, 



s(-l)l''l| 



<Ciexp(-C2^). 



(4.9) 



With this result we prove the following stronger version with strong concentration of Xpi^^ around 
q'^. Similar bounds on this quantity had previously been developed in the colouring model [6J to 
establish fast mixing of the block dynamics on tree and our proof is partially adapted from theirs. 

Lemma 4.2. When X > Ac, q^ < ^ andq'^q^ < l/(d— 1) there exist constants Ci(A,d), C2(A, d) > 
such that for s £ { + , — } and for large £, 

P (|Xp,,,, -q'\> exp(-Ci^)) < exp (- exp(C2^)) • 

Note that the condition ^^ < | is not necessary but simplifies the proof and holds in the regions 
of interest. 



Proof. We first observe that the Xp^i^s also satisfy a standard recursive relationship. If u G T^ and 
vi, . . . ,Vd-i are its children then the standard tree recursion for Gibbs measures of the hardcore 
model gives, 



Xv,, 



Anti(i-^^. 



(4.10) 



Now note that for any 5 > there exists £'{d, A, 6) such that for i > i' we have that X^/,s < q^ + 6 
since this is the case for even conditioning S^^i to be all or 1. Now for < L < i write Af/,/,s for 
the vector {X^^^^l^s : v G <S'p,l}. Observe that by recursively applying equation ()4.10p we can write 



X 



p,t,s 



9{^L,e,s)- 



Suppose that ?^l/,s,'^l e s ^^^ ^^° vectors which are equal except at some u G Sp^L- We will now 
estimate |5('^l,^,s) — g{X^ ^ ^)\. First consider one step of the recursion ()4.10p (i.e. the case L = 1). 
Let Ml, ... , Ufi-i be the children of p and suppose that u = ui. This implies that 



\giXi,£,s) - gi^i 



1 + Anti (1 - Xu.,i-i,s) 1 + Anti (1 - ^;/-i,j 



Now for a > we have that 



d 1 



dx 1+ax 



{1+ax) 



Y- If X > o then by a simple optimization we have 



that 'T .2 < |. Hence if ^ - L > i'{d, A, 1/15) then 



min{l - Xu,,e-i,s, 1 - <„^-i, J > 1 - (^7+ + j^) > 3 • 
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It follows that 



\9i^l,£,s) - 9iX[,i,s) 



i-x 



;.-!,. Ant2(i-^«./-iv 



< 



1-X„,,.i,, (1 + xAnt2 (1 - Xua-l,sW 

Y' I 



-dx 



u,l-l,s 3 

-dx 



^ \Xu,i-l,S 



Recursively applying this relation implies that for all L such that I — L > I' (d, A, 1/15), 



\9{Xlas) - 9iKi,s)\ < 



X 



u,£—L,s 



X,, 



ul—L,s\ 



By the Markov property of the configuration we have that the elements of XL,t,s are conditionally 
independent given ^^'^. Moreover, for u G Sp^i the Markov property also implies that Xu,i,s 

0) are strictly bounded away 



5.P Q^... TP/t^.P = 1) and 
from independent of L we have that by equation (j4.9p , 



depends on ^*'^ only through £_u'^ . Since 



u 



E 



IX, 



u,t,s 



s(-l)l"l 



e: 



i>p 



<C[exj>{-C2i). 



(4.11) 



Now choose some constant < r < 1 such that r\og{d — 1) — C2(l ~ r) < rlog(5/4) and set 
L = \ri\ . By Markov's inequality. 




-L,s 



A-^V 



> 



)' I er. 



< 



< 



< 



E exp E„, 



■J-&Sp L 



X 



ul—L,s 



A'^Y 






n^es,,^ K (exp ( 


exp((|)^) 
Xu,e-L,s - q'^-'^" 


] 1 ^^A 


exp((|)^) 
n.e5,,(l + eClexp(-C2(^-L))) 



exp((t)^) 



< exp (^{d - l)^eC; eM-C2{i - L)) - (^)^) 

< exp(-exp(C2^)) 



where the last inequality holds for large £ when < C2 < J' log |. Now if Xu/-l,s = 1 for all 

u G Sp^L then g{XL/,s) = 9* by the standard tree recursions. By equation (j4.1ip we have that if 

then 

,3 5 ,^ 



- ^4^ 



X, 



p,,s g1<0((---)^ 



and so the lemma holds taking Ci < '' log y| . ■ 

4.1. The measure on cry. For compactness of notation we will write the results of this section 
in terms of the plus phase but the analogous results will hold equally for the minus phase. Let Q^ 
denote measure on {0, 1} given by 

Q+(r/)=P^(ac; = r?|y(cT) = +). 

Lemma [33] shows that in expectation at least Q^(?y) behaves like Q^{r]). 



26 ALLAN SLY 

The graph G consists of G together with a collection of (d — l)-ary trees attached to U. Let 
Iri{G\ G) denote the independent sets on G\ G which are compatible with the boundary condition 
7]. Then the measure on the new part of G is given by 



^G^in\nM,u\Yi^) = + 



'(G\G)uu I -' v^y - T ; - — ^ u' - 1 • 

^a'e{o,iYG\G)uu^uy^u)^{cT'-(^r, iG\G)}'^ 



Now since Q'^{ri) « Q^{r]), at least in expectation, it will also be of interest to consider the measure 



Q^(-C/)1{.^^,,.^^(C\G)}A''^-^^' 



'(G\G)UUJ- , \a' \- 

^a'e{0,l}(G\G)uuWu{^u)^{cT'-Gl, iG\G)}'^ 

G\G u 

The graph {G\ G) Li U consists of {d — l)-ary trees of depth 2[^ log^-i n] rooted at each of the 
vertices of V and the leaves constitute U. For each t; G y let Ty denote the tree attached to v. 

Lemma 4.3. A configuration a G {0, l^^ ' distributed according to P* has the following 
properties: 

(1) The collection of projections {c^Tv^vev o,'''^ independent. 

(2) For each v G V^ the measure on ar^ is given by the projection of fi^ onto the first 
2[^ logf^^i n\ rows of the infinite {d — l)-ary tree. 

Proof. Since the trees T^ are disconnected and Q^ is a product measure the independence of the 
(TT„ is immediate. Verifying the distribution of o"t„ can easily be calculated via direct computation 
of the measure on the trees. However, we present a different proof which we feel better illustrates 
the replica method intuition underlying the result. 

Suppose that in the graph G we had that Q^{ri) = Q^{rj) holds exactly. Then the projection of 
the measure onto {G\G)DU is exactly given by P* . Note that this did not depend on the structure 
of G except through Qijirj). So consider the graph G* which consists on 2m' infinite {d — l)-ary 
trees whose roots we identify with U. Now take the Gibbs measure P^, on configurations on G* 
as follows: the measure is a product measure over the different trees and the measure restricted to 
an individual tree is p,~^ for trees rooted in U~^ and fi~ for trees rooted in U~ . 

Note that with this choice of graph and Gibbs measure the measure FQt{au G •) is Q^. Now 
construct G* by appending trees onto U identically as in the construction of G. The resulting 
graph is a collection of 2m disconnected {d— l)-ary trees rooted at the vertices of V. The resulting 
Gibbs measure ¥g* is a product measure over the trees of G* and restricted to individual trees 
it corresponds to the measure /*"'" for trees rooted in V+ and fi" for trees rooted in V~. By 
construction the measure on (G* \G*)UU is identical to P* which completes the proof. ■ 

Appending trees onto G to construct G has the effect of reweighting the projection of the measure 
on {0, 1}^. If we denote 

<Te/^(G\G) 
for f] G {0, l}*^ then 



¥g {^u = v\ y{(^) = +) 



E^G{o,i}^ QuivXv) 
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while 



{(^u = V) 



QUvMv) 



E^e{o,i}t^ QuivXv) 
Let B denote the set of configurations t] which have a large influence on ay by 



(4.12) 



B 



i; 



f] G {0, 1}^ : sup |PgK = l\au = r])- Q+K = 1)1 > exp(-2CiL^ log^.i n\) } , 



vev 



where Ci is as in Lemma 14.21 Observe that by the Markov property that the distribution of a^ 
depends only on rj through the leaves of T^ and hence that B is independent of G. By Lemma 14.21 
and a union bound we have that 



{au &B) < \V\ exp ( - exp ( 2(2 [- log^-i n\ 



We are now ready to prove Theorem 12.11 



(4.13) 



Proof. (Theorem 12. ip 

The number of vertices follows immediately from the construction. First choose 9 small enough 
so that 



a < mmj— — — -, 1. 



(4.14) 



51og(d-l)' 3 

We will first prove equation (|2.2p . Since ay is conditionally a product measure given ajj = r] we 
have that then for all -q' S {0, 1}^ that, 



Pg (o"y = V I 0-J7 = T/) _ ^ 



< 



-r-r Wg {(Ty = rj'^ \ au = 1]) ^ 
11 g+(a, = r?0 



vev 



Now ifr/G{0,l}'^\S then by the definition of B we have that 



^G 



[ay = 7?; I o-f/ = 7?) 



1 



<0(exp(-2CiL|log,_inJ) 



Qvi^v =Vv) 

It follows that if we take 9 according to (I4.14p then for large n and all r] G {0, 1}^\,B and rj' G {0, 1}^ 
then, 



Pg (o"y = r]' \ au = 1]) 



QtW) 



1 



< n 



-39 



since \V\ = 0{'nP). Hence we have that 

¥G{av = Tl'\Y = +) 



max 

v' 



QU^v) 



< n-^e ^ ^GicTu gB\Y = +) 



min^/ Q^(r/') 



(4.15) 



(4.16) 



Now since \V\ = 0{'nP) we have the inequality max.^/ {Qy{ri')\ = exp (0{'nP)\. To get a bound 
on Pg(o-{/ £B\Y = +) recall that 

T.r,eB^Gi^u = 'n\Y+)K{r,) 



G 



{au = v\ Y{a) = +) 



E^^e{o,i}^ ^G^^u = ri\Y = +) K{r]) 



Eve{o,i}^ ^Ci^vXv) 



(4.17) 
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Now by Theorem 13. 101 



while by Markov's inequahty 



{0,1}'^ '?e{o,i}c^ J 



(4.18) 



^Z+(7?H^)>iV^E^Z+(r?)K(r/)| ^ 0. 



k^eB 



r;GB 



Further, recall that by Lemma 13.31 for all r], 



EZt(r?) = (l + o(l))Q+(r?)EZ+ 



G 



G' 



(4.19) 



(4.20) 



and hence by equations (|4.12|) and (j4.13p 



E,sgEZt(r?)K(r?) _ ^^ ^ ^^^^^ E,eB IEQf^(^)'^(^) 



^»;e{o,i}t^^^(5(^)'^(^) 



(1 + «(!)) 



E,,6{o,i}f^IEQ+(r/)K(r/) 



^ 



= P* (at/ G i3) < |C/| exp ( - exp (2(2 [^ Iogrf_i n\ 

Combining equations KTf\\ . KT8\\ . ()i39|) . (Ii:20|) and (lOTTl it follows that for large n, 

Fciau eB\Y = +)<n\V\exp (- exp (^2(2 l| logd_i n 

Now provided that 9 satisfies (j4.14p then for large n, 

Fg{(Tu eB\Y = +)< exp(-n2^). 
Substituting this into (|4.16p 



(4.21) 



max 

r,' 



'G 



(av = v'\Y = +) 



1 



< 2n' 



-39 



QU^v) 

with room to spare for large n. The analogous statement for the minus phase 

PG(cTy = r/'|y = -) 



max 



Qvi'^v) 



1 



< 2n' 



-36* 



holds similarly and combining the two establishes equation (12. 2p . 
To establish (|2.ip we will show that with high probability 



^ r^ ^ G - 



^, 



G 



n Z, 



G 



n 



By equation ^^ EZt = (1 + o(l))EZ- and dMl) shows that for all r/ G {0, 1}^ we have that 
EZ^(?7) = (1 + o(l))EZ^(?7). In particular we have that since 



G 



and so 



EZ±= ^ EZ±(7?)^(r?) 
7je{o,i}'^ 

EZ+ = (l + o(l))EZg. 
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Now by Theorem 13. 101 



1 



[ Y. ZGin)<r,)<^¥. Yl Z^{rj)K{rj)]^0 (4.: 



Zt < —^EZ^ 1 = P I > Z^(v)k(v) < T-^E > Z^(v)k(v) I ^ (4.22) 



while by Markov's inequahty 



Z+ > -^/^EZ+ ) -^ 0. (4.23) 



1 
3 
Combining the previous two equations estabUshes equation (j2.ip and completes the proof. ■ 

5. Technical Condition 

Our last result is to verify Condition 11.21 in the case that A = l,d = 6. With ga, 13(^,6,8) = 
/(a,/3,7, (5, e) we use a computer assisted proof to show that ^^./S attains its unique maximum in 
the set ()1.6p at the point (7*, J*, e*) = (a^, /3^, q(1 — a— /3)) for (a, /3) in a neighbourhood of ip~ ,p'^)- 
These values are approximately p+ ^ 0.40831988,^" ^ 0.03546955 (see [9]). In [SHJ Lemma 6.3] it 
is shown that Qa^js is maximized as a function of e by taking 

1 r 



i{a, /3, 7, 5) = - [1 + a - /3 - 27 - v/(l - « - /3)2 + 4(a - 7)(/3 - <5)J ■ (5.1) 

It thus suffices to show that ga.i3{-, •) = 9a,(3i', ")£) is maximized at (a^,/3^) for 

0<7<a, 0<5<^. 
It is easy to establish that if 

r(a,/3,7,<5) = /i(a,7) + /'(/3,5) 

where 

f\a, 7) = H{a) + Hi{j, a) + i7i(a - 7, 1 - a) 

/2(/3, 5) = i/(/3) + Hii6,(3) + H^{p -6,1-13) 

then f < f* since / — /* is the log of a probability. Now as a function of 7, fi{a,j) is maximized 
at is maximized at 7 = q^. Similarly as a function of 6, /2(/5, (5) is maximized at (5 = /3^ and is 
increasing (resp. decreasing) in 6 for (5 < /?^ (resp. 6 > 13'^). Direct computation then shows that 

9p-,pA{P')Mp'^f) > 1-430 > 1.425 > fHp-,{p~f) + fip+, 0.015) 
9p',pA{P')Mp^f) > 1-430 > 1.414 > /i(p-,(p-)2) + /2(p+, 0.330) 

so for (a, /3) in a small enough neighbourhood of {p^ ,p^) and 6 G [0, igggl'-'CM)' /^l ^^-^"^ all < 7 < a 
we have that 

ga,l3{a,/3) > ga,f5{7,6) 

so it suffices to consider the set 

< 7 < a, -^ <6<^^ (5.2) 

- ' - ' 1000 - - 100 ^ ^ 

which we will denote T = T(a,/3). By |28l Lemma 5.1] the function ga,i3{-,-,i) has a stationary 
point at (a^, /3^) so the result will follow by showing that the Hessian matrix D'^ga^p{-, ■) is negative 
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definite in the region defined by ()5.2p . Tfiis in turn follows as we have that D'^g„^p+ is negative 
definite at ((p~)^, (p"'')^) (via a direct computation) and that 

detD^go,A-r)>0 (5.3) 



in the region defined by (|5.2p for {a, 13) sufficiently close to {p ,p~^). This is performed using a 
computer assisted proof. By [28^ Lemma 6.4] this determinant is given by 

,2. , ,, _ f df ^ da df \ f df ^ di df \ / df di df ^ ^ 



detZ) 5a,/3(T,^) 



+ 



+ 



i;;i)- 



\d'^j dj djde J \d'^S ' 86 d6de J \djd6 ' d"fd6de^ 
where the expressions for the partial derivative are given in |281 Lemma 6.2 and 6.4] as follows, 

df _ 6 6 5 

9^7 



2/3 + ,5 - 7 
6 



a — 7 
+ 



e "^ 1 - 2q + 



6 



df 
d'-fde 

K 



13-5- {a --f-e) 1-/3-7 
6 6 



+ 



a — 7 
6 



1 

7 



213 + 5 
6 



7-e 



a 



2/3 + 5-7 
6 



+ 



7- 
5 



+ 



/3-J-(a-7-e) l-p-^f 



13-5 -{a 



df df 



-7 
6 



+ 



2/3 + 5 

4 1 

p-5 ~ 5 

+ 



6 



and 



959e 975<5 1-2/3 + 5-7-e /3-(5-(a-7-e) 



/3-5 



5£__^ 

^7 ~ v'(l-"-/3)^ + 4(a-7)(/3-'^) 



05 



a 



7 



V(l-a-/3)2 + 4(a-7)(/3-5) 



Showing that the determinant is always positive is done using a computer assisted proof. Mathe- 
matica can perform interval arithmetic which given a function /i(-) and an interval [x, y] will return 
an interval containing the range of h([x,y]). This gives rigorous upper and lower bounds on the 
function including rounding in a conservative (i.e. rigorous) manner. This approach is slightly 
complicated by the fact that a couple of the terms go to infinity at the boundary of ()5.2p . We, 
therefore, do our estimates in a couple of stages. First let 



h^ 



df di df 
d'^5 d5 d5de 



Then 



max h-i(p ,p .■y,5,e) < max hi i p ,p' 

{■y,5)(iT 0<i<99,l<i<32 ' 



m a{i + 1) 
100^ 



100 



J_ J + 1 

100' 100 



< -17. 



where the last inequality is derived using mathematica's interval arithmetic. In particular we can 
always take hi to be negative for (a,/3) close to (p~,p^). 



We now analyse the term 2^— 



P-5 



6 



which appears in 



df I di df 



As 7 goes to a this term diverges so we need some further estimates on it before applying interval 
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arithmetic analysis. When < y < | we have that 



1 + ^y) <i + y 



and so since 



4(a-7)(/3-5) 



< 



Aa(3 



(l-a-/3)2 - (l-a-/3)2 
when {a, 13) is close to {p^,p^) we can take 



<0.19 



2 4(a-7)(/3-^) ^ / 4(a-7)(/3-^) 
5" (l-Q-/3)2 "V (l_a-/3)2 • 



Rearranging we conclude that 



V(l-a-/3)2 + 4(a-7)(/3-<5)- (1 - a - /3) < 



4(a-7)(/3-<5) 



5^(l-a-/3)2 + 4(a-7)(/3-<5) 
Now using this inequality and plugging in the definition of i we have that 



a — 7 — e = a — 7 

_ 1 r 
~ 2 



1 + a - /3 - 27 - V(l-a-/3)2 + 4(a-7)(/3-(5) 

V(l - a - /3)2 + 4(a - 7)(/3 - (5) - (1 - a - /3) 
(a-7)(/3-<5) 



5^(l-a-/3)2 + 4(a-7)(/3-5) 



and hence we have that 

4 



/3-(5 



6 



< 



-1 



a-7 ^(1 - a - /3)2 + 4(a - 7)(/3 - (5) a - 7 - e a-7 
for (a,/3) in a neighbourhood of {p~,p^)- It follows that 
df di df 



+ 



< 



6 



+ 



6 



+ 



6 



927 ^7 973e 1 — 2/3 + (5 — 7 — e 1 — 2a + 7 /3 — 5— (a — 7 — e) 1 — /3 — 7 — e 

1 1 



+ -1 + 



max{To^,a-7} max{^,7} 
/3-5 



V(l-a-/3)2 + 4(a-7)(/3-5)_ 
6 6 



+ 



6 



1- 2/3 + 5-7-6 /3-(5-(a-7-e) 1 - (3 - j - i 
where we denote the right hand side by ^(a, /3, 7, 6). Denote 



$(a,/3,7,(5) := 



xi,( R ew g/ I di df \ ( df di df Y 

vi'ia,p,7,oj ( ^2 J + 06d6deJ \d-fd6 d-f d6de J 



Since hi{p ,p'^ ,^,S,i) is negative throughout T we have that for (a,/3) in a small neighbourhood 
of (p~,]9+), 

max detD^^Q,^^ > max <I>(a, /3, 7, (5). 

(7,(5)eT ' (7,<5)eT 



32 ALLAN SLY 

Now applying a computer assisted proof using interval arithmetics we get that 



max ^(p ,p , 7, 5) > max <l> ( p ,p , 

(7,<5)eT 0<j<99,l<i<32 ' 



ai a{i + 1) 

Too' 100 



J_ j + l 

100' 100 



> 1500. 



By continuity of ^ this inequality also holds for (a,/3) in a small neighbourhood of {p ,p'^). This 
then establishes Condition 11.21 in the case that A = 1, d = 6. 
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